Detecting separable states via semidefinite programs 
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We introduce a new technique to detect separable states using semidefinite programs. This ap- 
proach provides a sufficient condition for separability of a state that is based on the existence of a 
certain local linear map applied to a known separable state. When a state is shown to be separable, 
a proof of this fact is provided in the form of an explicit convex decomposition of the state in terms 
of product states. All states in the interior of the set of separable states can be detected in this way, 
except maybe for a set of measure zero. Even though this technique is more suited for a numerical 
approach, a new analytical criterion for separability can also be derived. 



I. INTRODUCTION 

Entanglement is one of the most important resources 
for quantum information processing. It allows us to per- 
form such tasks as teleportation, secure quantum key dis- 
tribution, superdense coding and quantum computation 
to name a few [l| . Because of its central role a great deal 
of effort has been put into its characterization in the past 
few years. 

One of the most important questions, and also one very 
difficult to answer in general, is whether a given mixed 
state is entangled or separable. The decision problem 
associated with this characterization has been shown to 
be NP-hard [2(], so a simple practical procedure to an- 
swer this question is not likely to exist. To tackle the 
problem, many incomplete criteria have been developed 
in the form of either necessary or sufficient conditions 
for separability. Several of these criteria are based on 
relatively easy to verify properties that separable states 
must satisfy. If a given state fails such a test, it must 
be entangled. However, if the state passes the test, the 
result is inconclusive. This type of criteria includes the 
positive partial transpose (PPT) criterion d, 0| and its 
extensions!^, 0] , the reduction criterion Q and the range 
criterion [8[ among others. 

In this paper we introduce a new criterion for sepa- 
rability that works in a complementary way. If a given 
state passes a certain test, the state must be separable. 
Furthermore, a proof of this fact is given as an explicit 
convex decomposition of the state in terms of product 
states. The criterion is based on the fact that separable 
states preserve their separability property when a local 
map with certain properties acts on them. This local 
map need not be a physically implementable map, or 
even a positive map. The key point is that the search for 
such map can be implemented as a semidefinite program 
(SDP), which is a class of convex optimization problems 
for which efficient algorithms are known [9( . Semidefinite 
programs have found widespread application in quantum 
information [f| [l^, [HI Q3- Even though this new cri- 



terion is formulated in a way that is best suited for a 
numerical approach, it can also be used to derive a new 
analytical condition for separability. 

The paper is organized as follows. In Section I we 
present the basic idea of our approach by looking at the 
action of local maps on separable states. In Section II 
we introduce the concept of base states and discuss the 
properties they must have. Section III introduces an im- 
proved approach that can detect more separable states. 
Section IV deals with the problem of characterizing the 
set of states that are detected by this procedure. Sec- 
tion V shows how to extract a new analytical criterion 
for separability and serves as an illustration of the tech- 
nique. Finally, the conclusions are presented in Section 
VI. 



II. LOCAL MAPS ON SEPARABLE STATES 

Let po be a separable state on H.a'Si'Hb, where for now 
we will consider that dim Ha = dim7is = d. Then, by 
definition, we know that po can be written as a convex 
combination of product states, namely 
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with < pi < = 1 an d Pa* an d p^g are density 

matrices in TCa and Hb respectively. Suppose that we 
have a linear map A on density matrices over Hb, that 
satisfies A(p B ^) > for all i (where the inequality means 
the matrix is positive semidefinite (PSD)). Then, if we 
apply the map X ® A to po, the resulting state p is, up 
to some normalization constant, also a separable state, 
since 



p = 1 ® A(po) = X> Pa ® MPb) 



(2) 
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gives a decomposition of p as a convex combination of 
product states. Note that A need not be a positive map, 
since it is only required to be positive over a (typically) 
finite set of density matrices over TLb- Furthermore, 
since A may not even be normalized (in the sense that 
Tr[A(p^)] may not be equal to 1), after the appropriate 
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normalization the probabilities pi in the decomposition 
of p may be different from those in the decomposition of 
Po- 

Our approach to identifying separable states will be 
based precisely on this simple observation. Basically 
given a state a, if we find a map A satisfying the above 
mentioned criterion, and such that a = X ® A(po), we 
can assure that a is a separable state. Furthermore, we 
will be able to provide a proof of this fact in the form of 
an explicit decomposition of a as a convex combination 
of product states. We will see that the search for this 
map corresponds to solving a linear system and checking 
for the appropriate positivity requirements afterwards. 
Then we will show that we can slightly modify our basic 
observation ([2]) to enlarge the set of maps that would al- 
low us to detect separable states, and we will show that 
this new search can be cast as a semidefinite program. 

It is clear that for this approach to be useful, we need 
to start with a separable state po for which we know an 
explicit decomposition of the form (JTJ) . We will refer to 
such a state as a base state, since it is used as a "base" 
from which to reach other separable states. The useful- 
ness of this technique depends on how the state po is 
chosen, and we will start by discussing what properties 
it must have. 



III. BASE STATES 

Given a state a the question of whether it satisfies 
a = X ® A(po) for some map A can be answered by solv- 
ing a linear system of equations. Both a and po can be 
decomposed with respect to a basis for density matrices 
over Ha ® Hb as 

d 

PO = ^ ' Rmnrs 

\m)(n\ (g> \r)(s\, (3) 

m,n,r,s=l 

and 

d 

a= Smnrs\m){n\<g>\r){s\. (4) 

m,n,r,s — 1 

We can also define a basis on the linear space of maps over 
operators in Hb given by the set of maps {Mjki}fj f. i—i 
that satisfy 

A ijkl (\r)(s\) = \k)(l\6 ir 6 js . (5) 

The most general linear map can then be written as A = 

A ijk i , for some arbitrary coefficients 
Then the equation a — X® A(po) reduces to a system of 
linear equations given by 

d 

Smnkl ^ Rmnrs %klrs- (6) 
r,s — 1 

This is a system of d 4 equations with d A unknowns. It 
will always have a solution, but we also need to impose on 



this solution the positivity conditions A(p^ ) > 0,Vi. If 
the linear system © has a unique solution, we can check 
these conditions by computing eigenvalues. If the system 
is undetermined we can solve for some of the coefficients 
Xijki and write A = Ao +^2j Aj, where Ao is the fixed 
part of the solution and the sum groups the remaining 
terms. Thus the problem reduces to checking whether 
the map A is positive over the set {p^}- This can be 
written as a linear matrix inequality (LMI) , and so we can 
cast our problem in the form of the semidefinite program 
(SDP) 

min 
subject to 0, Aofpj?) + Ej xj ^j(Pb) > «■ (?) 

The size of the SDP will be determined by the dimension 
d and the number of terms in the decomposition |T]) of po . 
Caratheodory's theorem [l3| assures that any separable 
state can be decomposed as a convex combination of at 
most d +1 product states, so the size of the SDP remains 
polynomial in d. We will see later that we would have 
ample freedom to choose the base state po, and so we can 
keep the size of the SDP under control. 

The uniqueness of the solution of the system ([6]) de- 
pends an whether the matrix of the system is invertible 
or not. From ([5]) we can see that this matrix will de- 
pend only on the components R mnrs , and hence it is 
completely determined by the base state po. From this 
point of view it is a bit cumbersome to determine whether 
a given base state makes the system (j6|) invertible or not. 
Fortunately there is an alternative way of looking at this 
problem that connects with previous work and provides 
us with a nicer and more practical characterization of 
base states. 

First, let us note that for a fixed base state po we can 
regard the expression X® A(po) as a function that assigns 
to every map A the operator X ® A(p ). Note that we 
refer to it as an operator and not a state, since we have 
imposed no conditions on A, which may not be positive. 
If we choose po = \ijj + )(il> + \ (with = E»=i I") 

the maximally entangled state), then we know that the 
linear function /(A) = X ® A(\tfj + )(tfj + \) has very impor- 
tant and useful properties. It defines an isomorphism 
between the set of completely positive (CP) maps, map- 
ping states over Hb to states over Ha, and the set of 
states over Ha ®Hb- This is known as the Jamiolkowski 
isomorphism [Til ]. In (l5j D'Ariano and Lo Presti stud- 
ied whether other bipartite states besides the maximally 
entangled state will induce a one to one mapping be- 
tween CP maps and states. They found that there 
are indeed such states and derived a very simple test 
to identify them: a state po will induce such a map- 
ping if the operator po = (Epo) Tn E, is invertible, with 
E = J2ij I*J')0'*I the swap operator. They called such 
a state a faithful state. Furthermore, from this charac- 
terization it is easy to see that the set of faithful states 
is dense over the set of all states, and so we can have 
faithful states that are separable. However, we should 
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note that separable pure states cannot be faithful. In 
fact, if p = \x)(x\ (g) \y)(y\ (where \x) = Y,i x i\i) and 
\v) = TriVili}), then Po = \y)(x*\ ® \y*)(x\, which is 
clearly not invertible. 

If we use a faithful state as a base state, the function / 
is a one to one mapping between CP maps and positive 
scmidefinite operators. In terms of the linear system (J6j) 
this means that the system must be invertible since the 
relationship between states and maps must be one to one. 
Similarly, if the system ^ is invertible the state po must 
be faithful. It is clear then that if we use a faithful state 
as a base state, a larger set of separable states can be 
detected than if we use nonfaithful states. 



map 

er = [Kg, A B + A A (g>l] (po), (8) 

with Ab S Kb and £ Ka, where Ka — {Aa ■ 
Aa(Pa ) ) ^ 0,Vi}. The maps of the form J® A B +A A ®1 
also form a cone K A b that contains both Ka and Kb- 
If po is a faithful separable state, the map © sends all 
elements of K A b into the cone of separable states, and 
hence detects at least the same set of separable states 
than the map g alone. And finally, the search for a map 
in K A b that satisfies ((5J) for some state a is indeed a 
semidefinite program. 



IV. GEOMETRIC PICTURE AND IMPROVED 
ALGORITHM 

There is a nice geometric picture of our technique to 
detect separable states. First let us note that the set of all 
states is the intersection of the cone of positive semidef- 
inite matrices (known as the PSD cone) and the hyper- 
plane defined by the normalization condition Tr[p] = 1. 
The normalization does not affect the separability prop- 
erty of a state, so in order to make our discussion simpler 
we will sometimes analyze and prove results in terms of 
the cone structure of a set of states. It will be straight- 
forward to understand those results in terms of the ac- 
tual set of states satisfying the normalization condition 
Tr\p] = 1. 

Let us assume that we have a faithful base state po 
that is separable and has a convex decomposition given 
by The set of maps K B = {A : A(p { §) > 0,Vi}, 
is a cone (that is, a set closed under linear combinations 
with nonnegative coefficients) . The linear map g given by 
g(A) = X g)A(po), maps Kb into another cone contained 
in the set of separable states (which is itself a cone). It is 
not difficult to see that this image cone has nonzero mea- 
sure. First note that Kb contains the cone of CP maps 
which is in one to one correspondence with the PSD cone 
in Ha ®'Hb (the cone associated with bipartite density 
matrices) via the linear map / defined in the previous 
section. Both the cone of CP maps and the PSD cone 
are embedded in vector spaces of the same dimension d , 
and since the PSD cone has nonzero measure, it follows 
that so do the CP cone and Kb- Since the map g is linear 
and one to one, it will send a nonzero measure set into a 
nonzero measure set. 

We have seen that our approach reduces to solving a 
linear system and checking certain positivity properties 
of the solution. In particular if we use a faithful state 
as a base state, such linear system has a unique solu- 
tion, so no semidefinite program is required. However, 
we can improve our technique by considering a slightly 
more general type of mapping. Instead of considering 
g{A) = X® A(po) we can analyze the more general linear 



V. CHARACTERIZING THE SET OF 
DETECTED STATES 

The image of the cone Kb under the map g determines 
the set of separable states that can be detected for a par- 
ticular base state po- We have seen that the best choice 
for po is a faithful separable state, since in that case the 
mapping g is one to one and hence the image of Kb is as 
large as possible. This leads us to the obvious question of 
how big this cone is, or how much of the set of separable 
states can be detected in this way. So far we know that 
the set of detectable separable states (for a fixed faith- 
ful base state) is a cone of nonzero measure contained in 
the cone of separable states. We do not expect it to be 
the whole cone of separable states, since the separability 
problem is NP-hard and our characterization is polyno- 
mial. However we will show that we can detect in this 
way any state in the interior of the set of separable states 
S, except maybe for a set of measure zero. 

One set of states that we can assure will be detected 
is the set of separable faithful states in the interior 
of S, which is a set whose complement (with respect 
to S°) has zero measure. Let K + (p ) = {a : a — 
Xg)A(po), with A a positive map}. This is also a cone and 
it is contained in the cone of states characterized by the 
SDP discussed in the previous section. Even though this 
is a smaller set of states, it is enough to prove all the prop- 
erties we want to present. Note that it is clear that for 
any faithful state po (actually any state), po € K + (po), 
since we can just take A to be the identity map. But this 
observation is not at all helpful since we need to know if 
po is actually separable in order for this approach to be 
useful. Fortunately, for faithful states in 5° we can prove 
a stronger result: 

Theorem 1 Let a be a faithful state in S° . There is 
another faithful state po € S° such that a € (K + (po))°- 

Proof: Let a € S°. Consider the linear map h(A) = 
X <g> A(cr). By the continuity of this map, there is 
a neighborhood A of the identity map X such that if 
A e A => h(A) £ S°. Furthermore, since the iden- 
tity map is invertible, there is a nighborhood B C A 
such that all maps in B are also invertible. Let inv be 
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the inverse function, that sends any invertible map to its 
inverse. Then it is clear that inv — inv -1 . This together 
with the fact that inv is a continuous function tells us 
that inv send open sets into open sets. Consider then 
the open set B -1 — inv(B). Since I £ B and X -1 = X, 
then X 6 On the other hand, the identity map is 

on the boundary of the set of positive maps since an ar- 
bitrarily small perturbation of this map can send a PSD 
operator into a non-PSD operator. Therefore the inter- 
section of B -1 with the interior of the set of positive 
maps is a nonempty open set C. Now consider the open 
set C _1 = inv(C). Clearly, C _1 C A and every element 
A of C^ 1 is a map that satisfies (i) A 6 A and hence 
h(A) E S° , and (ii) A is invertible and A -1 is a positive 
map. Finally, consider the set h(C -1 ) C 5°. This is a 
nonempty open subset of S° , and so it must contain a 
faithful state po (different from a since the identity map 
is not in C _1 ), because faithful states are dense in the 
set of states. Then there is a map A -1 in C _1 such that 
pa = X <E> A -1 (cr). Therefore the map A £ C satisfies 
(j=I® A(/?o), with A positive. Since C is an open set 
contained in the interior of the set of positive maps, and 
the function 5(A) = X <£> A(po) is linear and one to one, 
a must belong to (K + (p ))°. □ 

The preceding theorem tells us that any separable 
faithful state can be detected using a different separa- 
ble faithful state as a base state. This comprises all of 
the set of separable states except for a set of measure 
zero (namely, the set of nonfaithful states and the states 
in the boundary of S that we have left out of this discus- 
sion.) We will see that some nonfaithful states are also 
detected by this procedure. To understand this we need 
the following lemma: 

Lemma 1 Let po be a faithful state in S. The state a 
'/, ® A(po) with A a positive map is also separable and 
faithful if and only if the map A is invertible. 

Proof: Consider the function /(A) = 1 ® A(cr). The 
state a is faithful if and only if this function is one to 
one. We can rewrite / as /(A) = I®(Ao A)(p )- If A 
is invertible then / is one to one, since /(Ai) = /(A2) 
implies that Ai o A = A2 o A because po is faithful, and so 
Ai = A2 because A is invertible. To prove the converse, 
assume that a is faithful and that A is not invertible, so 
that Ker(A) 7^ 0. We can choose two maps Ai and A2, 
Ai 7^ A2 such that they are distinct only over Ker(A). 
Then we have that Ai o A = A2 o A. Now we can write 

/(Ai) = T®Ai(o-) 

= I®(AioA)(po) 

= J®(A 2 oA)(p ) 

= X® A 2 (<r) =/(A 2 ), (9) 

which is a contradiction to the fact that / is one to one 
since a is faithful. Thus, A must be invertible. □ 

This result shows that there are nonfaithful states that 
are detected, since we can easily construct a noninvertible 
positive map A and hence the state X ® A(po) is not 



faithful and belongs to K + (po) for po faithful. However 
this is not enough to say that all nonfaithful states can 
be detected. Whether there are nonfaithful states that 
cannot be detected in this way is an interesting and open 
problem at this point. 

Another interesting result is given by the following 
lemma: 

Lemma 2 Let p\ and p 2 be faithful states such that P2 € 
K+{ P i). Then K+{p 2 ) QK+{ Pl ). 

Proof: Let a 6 K + {p 2 ). Then a = X (g) A(p 2 ) for some 
positive map A. Since p 2 € K + (pi), there is a positive 
map A' such that p 2 = X ® A'(pi). Combining the two 
equations we have that a = X ® (A o A')(pi). Since the 
composition of two positive maps is also a positive map, 
then g e K + {pi), and hence K + {p 2 ) QK + {p x ). □ 

An interesting picture of the set of separable states be- 
gins to emerge. By combining Theorem 1 and Lemma 2 
we can conclude that for any faithful state p 2 in S° we 
can always find another faithful state p\ also in S° such 
that K + (p 2 ) C K + (pi). We could then construct a se- 
quence of cones K + {p n ) such that K + (p n ) C K + (p n+ i). 
We do not expect this sequence of cones to contain all 
possible detectable states asymptotically, but rather to 
grow to some sort of maximal cone that will define a do- 
main inside the set of separable states. Since this domain 
cannot be the whole set S, we expect that there will be 
other domains constructed in the same way, using faith- 
ful base states that are not contained in any of the cones 
K + {p n )- These domains must cover the whole interior 
of S except maybe for a set of measure zero, comprised 
only of nonfaithful states. 

To summarize, the procedure for detecting separable 
states could be applied as follows. Given a state a we 
want to check for separability, we first generate a random 
separable state by randomly generating the terms in the 
decomposition 

i 

i.e., randomly choosing the number of terms in the de- 
composition, the probabilities pi and the local density 
matrices p^ and p^ . We can simplify this step by not- 
ing any separable state can be written as a convex combi- 
nation of pure product states with at most d A + 1 terms. 
Then we check whether this random separable state is 
faithful or not by checking if the operator p = (Ep ) TB E 
is invertible (with E = £\. \ij)(ji\ the swap operator.) 
Since the set of nonfaithful states has zero measure, our 
random state will be most likely faithful. Then, we can 
just solve the enhanced SDP discussed in Section HVl and 
if a feasible solution exists, we have proven that a is sep- 
arable, and can furthermore construct an explicit convex 
decomposition of it in terms of product states. If the 
SDP turns out to be unfeasible, it means that a could 
be entangled, or that we need a different faithful base 
state. We could then select a new faithful base state and 
proceed as before. 
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Since each faithful state can help us to identify many 
separable states, it could be useful to save them for use 
with other states. We can then build up a table of faithful 
separable states for each dimensionality, so we can bypass 
the random generation step in the future. Furthermore, 
this table can be optimized by discarding any state p2 
that belongs to K + {pi) for some pi also in the table, since 
according to Lemma 2 it will not detect any new states. 
Another interesting and open problem is whether we can 
make this table finite and still detect most of the states 
in S° except maybe for a set of arbitrarily small (but 
not zero) measure. That finite number of base states will 
probably depend on the dimensions of the problem (and 
the measure of the undetected set), but this approach 
might turn out to be a useful tool for low dimensional 
problems. 

A. Extension to <1a 7^ cis 

Our approach to detect separable states can be ex- 
tended to the case in which the dimensions of the spaces 
Ha and Hb are not equal, with a few minor modifi- 
cations. To see how this works we first need to state 
the Jamiolkowski isomorphism in its more general form. 
Let Ha and Hb be two spaces of dimensions cIa and 
ds respectively. We will denote by B(H) the set of 
bounded, PSD operators over the Hilbert space H. Let 
C(B{H A ),B(H B )) be the set of CP maps from B{H A ) 
to B{H B ). Then the function F : C(B(H A ),B(H B )) -» 
B(H A ®H B ), defined by 

F(A) =1®A(\iP + ) aaaa (^ + \), (11) 

where \ip + ) AA is the maximally entangled state in 
Ha <& Ha, is an isomorphism between CP maps in 
C{B{H A ),B{H B )) and PSD operatos in H A ® H B . 

Now let us define another function F : 
C(B(H A ),B(H B )) -> B(H A ®H B ), given by 

F(A) =1® AO*), (12) 

with po a faithful separable state in H A (g> Ha- It is clear 
that F(A) is also a separable state in Ha ®Hb provided 
that A satisfies the same positivity conditions we required 
in the case cIa = ds- By the same isomorphism, applied 
to states in Ha <8> Ha, we know that there must be a CP 
map A : B(Ha) — » B(Ha) such that 

Po = l®A (\<ip + ) AAAA {ip + \). (13) 

Furthermore, following the same reasoning as in the proof 
of Lemma 1, we conclude that the map Ao must be in- 
vertible (note that this inverse will not be in general a 
CP map.) And using the same argument once again, it 
is not difficult to show that the function F must be one 
to one. By the same arguments of the case cLa = ds, we 
can then conclude that the set of detected states is also 
of nonzero measure when dA 7^ ds ■ The only difference 



between these two cases is that in the latter the faithful 
states must be chosen from a space different than the one 
in which the states to be detected live. 



B. Multipartite states 

It is interesting to see what happens with this ap- 
proach when we apply it to multipartite states. Even 
though the basic property that allows us to connect two 
separable states is still valid, the results are not as use- 
ful because the set of detected states has zero measure 
even for faithful states. Let us briefly show the reason 
in the tripartite case. Let po be a separables state in 
Ha <8 Ha d> T~t A , with dimH A = d. Clearly, any state of 
the form a — [T (8 1 ® A](po), is separable if A satisfies 
the required positivity constraints. The cone of separable 
tripartite states is embedded in a vector space of dimen- 
sion d 8 . But the set of linear maps A is embedded in a 
space of dimension d , so even if po is a faithful state, 
the mapping of A to a produces a set of dimension d 4 
embedded in a space of dimension d 8 . Hence, the image 
of this mapping (which is the set of detected states) has 
zero measure with respect to the cone of tripartite sepa- 
rable states, and this makes this approach not useful in 
practice. One solution to this problem will be to consider 
transformations of the form a — [X <£> Aa <8> As](po) but 
then checking for the existence of the maps A^ and A# 
becomes a nonlinear problem that cannot be solved using 
SDP. 



VI. ANALYTICAL CRITERION 

This technique can also help us to derive some analyt- 
ical results about separability. By choosing a particular 
base state, and using some results on the characterization 
of positive maps [f| , we can obtain a sufficient condition 
for separability in terms of the eigenvalues of an associ- 
ated state. 

Consider the following family of states 

p(A) = (l-A)/ A B + A|^ + )(^ + |, (14) 

where < A < 1, Iab is the maximally mixed state in 
Ha ® Hb (i.e., the identity matrix times 33, since we 
are considering dimHA = dim7is = d), and \ip + ) = 
^Eil^)^B is the maximally entangled state. These 
states are referred to as the isotropic states (since they 
are invariant under U ® U* transformations) and are 
known to be separable for < A < ^L-, and entangled 

otherwise. We will choose po = p{j^i)- This state is in 
the boundary of S and it is not difficult to check that it 
is faithful. 

Let a be a state that satisfies 

a=l®A{ Po ), (15) 
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for some positive map A. This is equivalent to 

a=l(g>A(^I AB )+l®A(j^\^+)(ib+\). (16) 

Let us impose a restriction on A and require A (7b) = Ib- 
We will see later that this restriction can be eliminated, 
but for now it simplifies the reasoning. Then we can 
rewrite ([TBI) as 



o = j^jIab + A(|^+)(V, + |). (17) 

The restriction on A means that, for consistency, we must 
have Tr^u] = Ib- Again, this restriction can be lifted 
and we will see that later on. We now write 



-I 



d+1' 



AB 



:Z®A(|V + }(^ + |), 



where we have rescaled the positive map, so A = ^pjA. 

From the Jamiolkowski isomorphism we also know that 
any operator Z of the form 



Z=X® A(|^+)(^+|) 



(19) 



with A a positive map, must be positive over the set 
of separable states. This is equivalent to the associated 
bihcrmitian form 



E z = (xy\Z\xy) 

ijkl 



(20) 



being nonnegative. In j| a sequence of tests was intro- 
duced (each test implementable as a SDP) that could 
prove the positivity of the form Ez, and hence the exis- 
tence of a positive map satisfying (|19p. 

Consider our state a that satifi.es Tr^ [a] = Ib and let 
us assume that we can prove that the bihermitian form 
associated with the operator (a — -jt^Iab) is positive, 
using the techniques from || . Then, by the Jamiolkowski 
isomorphism, the map A that satisfies 



d+1 



I A B=l®A(\i> + )(ip + \), 



(21) 



is positive. By tracing over A on both sides it is not 
difficult to show that A(I B ) = j^[Ib, so we can write 



d+1 



'AB 



d+1 



J®A(|^+)(^ + |), (22) 



where now A satisfies A(I B ) = Ib- Since Iab = I A® Ib 
we can write Iab =1 ® A(Iab) ancl then we have that 

o = 1®A{-^—Iab)+I®K^-t\^ + )^ + \) 
d+1 d+1 

= 1®A{ Pq ), (23) 



Assume now that we have a state a that satisfies 
Tr^[cr] = <j B > 0, so that a B is invertible. Note that 
any state in the interior of S will have this property. We 
can construct the state 



1 



(l A (g>a B 2 )a(l A 



(24) 



that satisfies the condition Tr^fa] = I B (with 1a the 
identity matrix in d dimensions). It is clear that if a is 
separable, so is a. 

According to our discussion above, we can prove sep- 
arability of a by proving positivity of the bihermitian 
form associated with the operator a — -A^Iab- Clearly, 
a sufficient (but not necessary) condition for this to be 
true is the that the operator a — -A^Iab itself be pos- 
itive semidefinite. If Xi are the eigenvalues of a, this is 



(18) equivalent to requiring that min^ > 



have the following corollary: 



3(3+1) 



We then 



Corollary 1 Let a be a state in Ha®'Hb with dim?^ 

dim7Ys = d, that satisfies Tr^u] = a B > 0. Let a 
_ i _i 

^(1a ® cr B 2 ) cr (1a <8 o B 2 ), with eigenvalues Xi- If 



min, A,- > 



1 



d{d+l) 



(25) 



then a is separable. 



Furthermore, for any state a shown to be separable as 
a consequence of Corollary 1, we can provide an explicit 
convex decomposition in terms of product states. From 
the Jamiolkowski isomorphism we can also extract an 
explicit expression of the map A in terms of the hermitian 
operator Z that satisfies Z = X ® A(|^ + )(t/> + |). In our 
case, given an orthonormal product basis {|u)ab}, we 
can write 



Hp) 



E<*j|(a 

ijkl 



1 



lAB)\kl)(i\p\k)\j)(l\, (26) 



where a is given by (j2"4"|) . Now we just need a decompo- 
sition of our base state 



Po 



d+1 



-I 



1 



AB 



d+1 



(27) 



In [16| one such decomposition of po w as given and we 
will follow that construction. First, let us define a vector 
z = (zi, . . . , Zd), whose components Zj take on the values 
±1 and +i. To each vector z we associate a pure state 



I*. 



d 



(28) 



with {\j)} the canonical basis. There are 4 d such vectors, 
and hence 4 d states |$ z ), although only 4 d_1 are distinct 
in that they differ by more than a global phase. Now we 
define a product state in Ha ®"Hb given by 



which shows that a is separable. 



p z = |*.)<*.|®|fc..><*..|. 



(29) 
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The ensemble consisting of all 4 d of these states with the 
same probability gives the density operator 

pE^ z = 4^ E w2 z i z t z i z m \ \j){k\®\l){m\. 

z jkhn \ z / 



(30) 



Since 



2J Z zlZiZ m = 4 d (5j k Sl m + SjlS km - SjkSl m Sjl), (31) 



it follows that 



-^$>=/Ai 3 +~iv + ><v + i - ^E u>oi ® u)oi 



rf 2 



(32) 



Multiplying by and rearranging terms we get 

p« = ^hT,p* + j^ d hm®\m, 03) 

7 = 1 



and hence we have 
d 2 1 



,/ + rid E ® ^a(|$ b .)(^- i vi + 

z 

+7^ryEii)oi®4A(U)oivi (34) 

7 = 1 



where <jb = Tr^c] and A is given by (|26|) . Since the 
map A is positive, Eq. (l34|) gives an explicit convex de- 
composition of a in terms of product states, certifying its 
separability. 

It is interesting to compare this corollary to the well- 
known result of Gurvits and Barnum [l7T ] that charac- 
terizes the biggest ball of separable states centered on 
the maximally mixed state. Their result says that if 
1 1 cr — I AB |||< gg7gg3Y) ' then the state a is separable. 
Now consider a state that satisfies Tr^[cr(e)] = Ib in the 
basis in which it is diagonal, and that 



cr(e) = diag(e + 



1 



d(d + l) 



A,..., A), (35) 



where A = d ^+i) + ^ an d both e and 5 are positive. By 
normalization we must have d ^ +1 ^ + e + (d 2 — 1)5 = 1, 
and hence 



S=(l-e 



d + V d 2 - 1 ' 



d+l ■ 



(36) 



If these 



Since we need 5 > 0, we need e < 1 
conditions are satisfied, the state a(e) is separable. To 
see what the Gurvits-Barnum criterion says about such 
a state, we define the function 



Me) =|l vV-Iab \\l 



1 



d 2 {d 2 - 1)' 



(37) 



If f(d, e) > for < e < 1 - ^j, then the state cr(e) is 
not shown as separable by the Gurvits-Barnum criterion, 
but it is by ours. It is easy to check that 



M,i 



d-2 



d + l' d 3 -d 



> 0, yd > 3. (38) 



Since f(d, e) is a continuous function of e, this shows that 
for d > 3, there is always a range of values of e for which 
the state cr(e) is shown to be separable by our technique 
but not by the Gurvits-Barnum criterion. 



VII. CONCLUSIONS 

In this paper we have introduced a new technique to 
detect separable states. The idea is to show that two 
states, the first of which is known to be separable and 
has certain properties, are connected by a local map that 
preserves separability. This local map is required to be 
positive on a set of local states determined by the convex 
decomposition of the first state in terms of product states. 
If a map with the required properties is found connecting 
the two states, the second state is proven to be separa- 
ble. The key point is that searching for this connecting 
map reduces to a semidefinite program (SDP) which can 
be implemented efficiently. Furthermore, if the second 
state is proven to be separable, the solution of the SDP 
provides the required map and an explicit proof of sepa- 
rability of the second state can be constructed in the form 
of a convex decomposition in terms of product states. 

It is important to note that the local map in question 
needs not be positive over all states. This is an advan- 
tage, since characterizing the set of positive maps is an 
NP-hard problem, while the wider set of maps we are 
interested in (which includes the set of positive maps) is 
easier to characterize and allows us to reduce the prob- 
lem to a SDP. However, restricting ourselves to the set 
of positive maps makes it easier to prove and understand 
some properties of this technique. In particular, coupled 
with some previous work aimed at characterizing positive 
maps via an infinite hierarchy of conditions [5j , it allowed 
us to extract a new criterion for separability based on the 
spectral properties of an associated state. 

We have shown that all states in the interior of the 
set of separable states can be detected by this technique, 
except maybe for a set of measure zero. Our approach 
proceeds by choosing a random separable state (with a 
certain property known as faithfulnes [HI) that we call 
base states, and implementing a SDP to check whether a 
given state is separable or not. The set of states detected 
using a particular base state is a convex subset of the set 
of separable states. By choosing different base states we 
are able to detect more and more separable states. We 
could keep a table of these states for future use much in 
the same way we keep the information about entangle- 
ment witnesses. Even though an infinite number of these 
states are needed to characterize all separable states, it 
may well be the case that a finite number will be enough 
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to characterize almost the whole set, except maybe for 
a set of very small (but not zero) measure. This is still 
an open question. Also open for further research is the 
structure of the set of separable states that cannot be de- 
tected using this technique. This structure seems to be 
related to the structure of the set of non invertible maps. 

Finally, it is worth stressing that, as many other sep- 
arability criteria, our technique is not complete, and 
should be complemented with other ways of analyzing 
entanglement and separability. One feature that should 
be pointed out is that our criteria either proves separa- 
bility of the state (and provides a certificate of it in the 
form of an explicit decomposition) or fails. Other criteria 
usually show that the state is entangled or fail. These cri- 
teria generate an approximate characterization of the set 
of separable states from the outside. Our approach works 
from the inside. In this respect it is similar to the algo- 



rithm presented in [l8| , that generates a sequence of con- 
vex sets included in the set of separable states for which 
membership can be easily checked. These two types of 
criteria (inside and outside characterizations) should be 
used together to increase the chance of correctly identi- 
fying the entanglement properties of a given state. 
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